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SIMPLIFIED  METHOD  FOR  CALCULATION  OF  COMPRESSIBLE 
LAMINAR  BOUNDARY  LAYER  WITH  ARBITRARY 
FREE- STREAM  PRESSURE  GRADIENT 
By  George  M.  Low 


SUMMARY 

The  Karman-Pohl hausen  integral  method,  as  applied  to  compressible 
laminar  boundary  layers,  was  simplified  by  an  analysis  similar  to  the 
incompres sib le- Holstein- Bohlen  method.  Although  this  simplification  is 
helpful  for  many  calculations,  it  is  of  greatest  value  when  conditions 
at  the  edge  of  the  boundary  layer  are  known  from  experimental  measure¬ 
ments.  The  analysis  was  conducted  under  the  assumptions  of  a  Prandtl 
number  of  1,  zero  heat  transfer,  and  a  linear  viscosity- temperature 
relation;  velocity  profiles  are  approximated  by  a  fourth-degree 
polynomial . 

Results  are  presented  so  that  velocity  and  temperature  profiles, 
momentum  and  displacement  thicknesses,  and  wall  shear  stress  can  be 
calculated  for  flows  over  two-dimensional  bodies  with  arbitrary  free- 
stream  velocity  distributions.  The  results  are  also  applicable  to 
flows  over  three-dimensional  bodies  with  axial  symmetry  through  the  use 
of  Mangier' s  transformation. 


INTRODUCTION 

Present-day  theory  of  compressible  laminar  boundary  layers  permits 
the  aerodynamicist  to  make  satisfactory  boundary- layer  calculations  for 
flows  with  zero  pressure  gradients.  Exact  solutions  of  the  boundary- 
layer  equations  for  flows  with  streamwise  pressure  gradients,  however, 
exist  only  for  several  specific  pressure  distributions  (reference  l). 
Solutions  of  the  laminar-boundary- layer  equations  for  flows  with  arbi¬ 
trary  pressure  gradients  are  of  ever  increasing  importance  because  it 
has  been  suggested  that  at  sufficiently  high  altitudes  and  speeds  lami¬ 
nar  boundary  layers  can  exist  at  very  high  Reynolds  numbers  (refer¬ 
ence  2),  and  because  the  pressure  gradients  encountered  over  most  wings 
and  bodies  are  not  approximated  by  any  of  the  aforementioned  exact 
solutions. 


2 


NACA  TN  2531 


The  lack  of  exact  solutions  for  flows  with  arbitrary  streamwise 
pressure  gradients  necessitates  the  use  of  approximate  solutions.  The 
KarmAn-Pohlhausen  method  has  long  been  found  to  yield  satisfactory  solu¬ 
tions  for  incompressible  flows,  provided  that  the  separation  point  is 
not  approached.  Howarth  (reference  l)  has  recently  applied  this  method 
to  compressible  flows. 

The  solution  of  the  incompressible  or  compressible  boundary- layer 
equations  by  the  Karman-Pohlhausen  method  is  complex  because  it  involves 
the  numerical  solution  of  a  cumbersome  differential  equation,  and 
because  it  involves  second  derivatives  of  the  free- stream  velocity, 
which  cannot  be  determined  with  any  accuracy  if  the  free -stream  velocity 
is  known  only  from  experimental  measurements.  For  incompressible 
fluids,  these  difficulties  have  been  overcome  by  the  Holstein-Bohlen 
simplif i cation  (reference  3),  which  eliminates  the  second  derivatives 
of  the  free-stream  velocity  and  simplifies  the  differential  equation. 

In  the  present  report  of  an  analysis  made  at  the  MCA  Lewis  labo¬ 
ratory,  the  Holstein-Bohlen  simplification  is  applied  to  the  Karman- 
Pohlhausen  method  for  compressible  fluids.  The  momentum- integral  equa¬ 
tion  is  derived  using  the  transformation  of  Howarth  (reference  l) 
together  with  a  more  general  viscosity- temperature  relation  than  that 
used  in  reference  1.  The  greater  part  of  the  analysis  presented  fol¬ 
lows  along  the  lines  of  reference  1,  but  the  emphasis  herein  is  placed 
on  simplifying  the  final  results  of  the  analysis  so  that  they  can  easily 
be  applied  to  any  practical  problem. 

Although  the  present  analysis  is  carried  out  for  two-dimensional 
flows,  it  can  be  applied  to  three-dimensional  flows  with  axial  symmetry 
by  using  the  transformation  of  Mangier  (see  appendix  A).  Under  this 
transformation,  a  three-dimensional  body  with  a  given  pressure  gradient 
is  transformed  to  a  related  two- dimensional  body.  The  results  of  the 
present  study  can  be  applied  to  this  two-dimensional  body. 


SYMBOLS 

■  following  symbols  are  used  in  this  report: 
speed  of  sound 

proportionality  factor  used  in  equation  (l) 
specific  heat  at  constant  pressure 
function  defined  by  equation  (28) 


The 


f(Mx) 
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fp(X)  function  defined  by  equation  (28) 

fg(X)  function  defined  by  equation  (28) 

f(rj)  function  defined  by  equation  (20) 

h  increment  of  lengtl^ 

i  summation  index 

K  arbitrary  constant 

k  thermal  conductivity 

M  Mach  number 

n  transformation  variable  defined  by  equation  (9) 

Pr  Prandtl  number,  pcp/k 

p  static  pressure 

r  radial  coordinate 

S  Sutherland^  constant  (216°  F) 

t  static  temperature 

u,v  velocity  in  x-  and  y-directions,  respectively 

x,y  cartesian  coordinates  measured  along  body  and  perpendicular  to 
body,  respectively 

Z  function  defined  by  equation  (25) 

X  ratio  of  specific  heats 

A  boundary-layer  thickness  in  x,n  coordinate  system 

A*  function  defined  by  equation  (19) 

5  boundary- layer  thickness  in  x,y  coordinate  system 

5*  boundary- layer  displacement  thickness 
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©  function  defined  by  equation  (l8) 

9  .  boundary- layer  momentum  thickness 

X  Pohlhausen  parameter 

H  absolute  viscosity 

v  kinematic  viscosity 

p  mass  density 

x  shear  stress 

0  transformed  stream  function 

t|t  stream  function 

Subscripts: 

0  conditions  at  stagnation  point 

* 

1  conditions  at  outer  edge  of  boundary  layer 

s  standard  condition  * 

w  condition  at  solid  boundary 

ASSUMPTIONS 

The  following  assumptions  are  made  in  addition  to  the  usual 

boundary- layer  assumptions  in  order  to  make  the  problem  amenable  to 
solution: 

(1)  The  Prandtl  number  of  the  fluid  is  equal  to  1.  It  has  been 
shown  in  reference  4  that  this  assumption  has  no  appreciable  effect  on 
flat-plate  skin- friction  calculations  if  the  fluid  is  air  (for  M<  5), 
but  the  assumption  does  limit  the  accuracy  of  heat- transfer  calcula¬ 
tions  and  temperature  distributions. 

(2)  The  heat  transfer  at  the  solid  boundary  is  equal  to  zero. 

This  assumption  should  be  reasonable  for  many  practical  applications. 
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(3)  The  viscosity  and  temperature  are  related  linearly  hy  the  fol¬ 
lowing  expression: 


Chapman  and  Rube  sin  (reference  5)  have  shown  that  solutions  of  the 
boundary-layer  equations  based  on  equation  (l)  agree  well  with  reality 
for  flat-plate  flows  if  the  constant  C  is  determined  by  matching 
equation  (l)  with  Sutherland^  relation  at  the  solid  boundary  so  that 


A/V(ts  +  s) 

v  (y+  s) 


(2) 


This  assumption  should  also  be  reasonable  for  the  case  of  flows  with 
streamwise  pressure  gradients  because  the  wall  temperature  t^  is  con¬ 
stant  as  a  consequence  of  assumptions  (l)  and  (2). 


ANALYSIS 

Compressible  laminar-boundary- layer  equations.  -  The  momentum  equa¬ 
tion  describing  the  flow  in  the  boundary  layer  is 

„  Su  c>u  _  1  dp  i  S  f  du\  ,  s 

U  3x  +  cty  ~  "  p  dx  * :  p  ciy  V  (3) 

The  equation  of  continuity  is 

(pu)  +  ^  (pv)  =  0  (4) 

The  flow  of  heat  is  described  by  the  energy  equation,  which  is 

pcp(us  +  TH)-ui  +  |;(ltlt)  +  ^(lj)2  <5> 

A  particular  solution  of  the  energy  equation  can  be  obtained  by  multi¬ 
plying  equation  (3)  by  pu  and  adding  the  result  to  equation  (5),  and 
by  letting  the  Prandtl  number  equal  1.  This  solution  is 


2 

and  corresponds  to  the  case 
ary  is  equal  to  zero. 


+  Cp  t  .=  constant  (6) 

where  the  heat  transfer  at  the  solid  bound- 
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In  terms  of  the  transformed  coordinate  system,  a  function  $(x,n)  can 
be  defined  as 


^  2 

t(x,y )=V^"(^)  ^x,n' 


From  equations  (id),  (ll),  and  (12), 


ps  L  /PC  dn 


Ps  V  dx  dn 


s) +  s'  )_ 


and,  as  shown  in  reference  1,  from  equations  (l),  (7),  (8),  (10), 
and  (13),  the  momentum  equation  becomes 


on  Sxdn  ”  ox  ^n2 


du 
U;1-  dx 


Jr  Vs 


The  temperat ure  term  in  this  equation  can  he  expressed  in  the  form 


from  which 


M  B2gi  & 


n  dndx  dx  ^n2 


2a!2  5n2 


1- 


S30 

S  ir? 


With  the  exception  of  the  term  within  the  bracket,  equation  (16)  is 
identical  in  form  to  the  incompressible  momentum  equation. 

Momentum- integral  equation.  -  The  momentum- integral  equation  is 
obtained  by  integrating  each  term  of  equation  (16)  from  the  solid 
boundary  to  the  outer  edge  of  the  boundary  layer  (see  reference  l) : 


2  d© 

ui  es +  u 


13F  8  (2‘¥)  -  (t  +  ^  «i2)  *  (i)„ 
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where 


The  auxiliary  functions  0  and  A*  are  similar  to  the  momentum  thick¬ 
ness  6  and  the  displacement  thickness  5*  in  the  physical  coordinate 
system. 


Modified  Karman-Fohlhausen  solution.  -  In  order  to  solve  equa¬ 
tion  (17)  "by  the  Kdrman-Pohlhausen  method  it  is  necessary  to  introduce 
the  variable  T)  =  n/A  and  to  assume  that  u/u^  is  a  polynominal  in  tj. 

As  in  the  method  of  Pohlhausen,  this  polynomial  is  taken  to  he  of  the 
fourth  degree.  The  coefficients  of  the  polynomial  are  determined  from 
the  following  boundary  conditions: 

n  =  0:  u  =  0 

-ui  IT  (x  +  ¥  Mi2)  -  vs 

n  =  A:  u  =  up 

Su  a2u 

^n  Sn2 

Thus  the  following  form  is  obtained  for  the  velocity  profile: 

^  =  f(il)  =  2TJ--  2t)3  +  T]4  +x[j  T)  (l-Ti)^J  (20) 

where  the  parameter  X  is  defined  as 

<*> 


d30  _  ^2u 

8n3  s 


a 
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Transformation  to  the  ^-coordinate  and  subsequent  integration  in  equa¬ 
tions  (18)  and  (19)  yield 

«.  ..  ; 

®=3b(57  -I'll)  <22> 


A*  =  A 


(56  - 1  \ 

V  120  ) 


and  it  is  found  that 


(*A  q/,  ii 

Uj„  =  a  (2  *  e) 


Equations  (22),  (^3),  and  (24)  could  be  applied  to  equation  (l7) 
directly  and  the  result  would  be  a  differential  equation  for  X  (ref¬ 
erence  l).  A  simpler  solution  can  be  obtained  by  applying  the  substi¬ 
tution  of  Holstein-Bohlen: 

Z  =  02/v  s  (25 

or,  from  equations  (22)  and  (2l), 


,7  X  5X_ 
»  ~_3  144, 

I+TM!2) 


Equation  (l7)  is  multiplied  by  ©/u-^,.  and  equations  (25)  and  (26)  are 
substituted  in  the  resulting  expression  to  yield 


ll  dZ 


f  S  +  *iW l_-  *<«i>  tjJ-t 

Then,  from  equations  (22),  (23),  and  (24), 


l _ 0_  /du\ 

J  ~  ui  vWv 


-  ^  £fl(X)  f(Ml)  +  f2^)J 
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•where 


and 


(28) 


(29) 


Equations  (27)  and  (29)  can  now  he  integrated  simultaneously,  and  X 
can  be  obtained  as  a  function  of  x  by  a  method  of  numerical  integra¬ 
tion  given  in  detail  in  appendix  B.  The  functions  f^(x)  and  fg(\) 
are  graphically  presented  in  figure  1  and  are  tabulated  in  table  I; 
f(Mi)  is  tabulated  in  table  II.  Boundary- layer  variables  are  still  in 

the  transformed  system  of  coordinates^  however 9  and  a  transformation  to 
the  physical  (x,y)  coordinates  is  required. 


Inverse  transformation  to  physical  coordinates.  -  The  distance 
normal  to  the  surface  can  be  obtained  from  equation  (9): 


or 


A 


m 
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where  A  is  obtained  from  equation  (2l).  The  boundary- layer  thickness 
is  determined  by  letting  t|  equal  1: 


+  Mx2  (-0.0001X2  -  0.0094X  +  0.4175) 


— 

) 


The  momentum  thickness  is  defined  as 


(33) 


and  the  displacement  thickness  is 


The  vail  shear  stress  tw  is  defined  as 
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The  separation  point  is  defined  as  the  point  where  the  shear  stress 
vanishes,  or,  from  equation  (36),  where  X  =  -12.  It  is  emphasized, 
however,  that  the  Karman-Pohlhausen  solution  has  been  found  to  be  quite 
inaccurate  in  the  neighborhood  of  a  separation  point. 


SUMMARY  OF  METHOD 

The  results  of  this  study  will  be  summarized  in  the  order  they  are 
needed  for  the  solution  of  a  particular  problem.  It  is  assumed  that  all 


free- stream  conditions 


( 


du- 


ul>  dx  ’  Ml>  - t 


1>  Pi) 


are  known  functions  of  x. 


The  standard  conditions  (vg,  pg,  tg)  are  selected  as  mean  values  of  the 

free- stream  conditions.  Next  X  is  determined  as  a  function  of  x  by 
solving  equations  (27)  and  (29)  numerically: 


where 


dZ 

dx 


=  JL  [fi(x)  f (%)  +  f2(x)] 

fi(x) 


dui 

dx 


( '^  'V) 


f2(X>  $  5E  (3?  -  fs)(2  -  1 


f(Ml)  = 


Mi  ‘ 


2  +  (r-1)  mx2 


(27) 


(29) 


(28) 


A  suggested  procedure  for  solving  equation  (27)  is  presented  in  appen¬ 
dix  B.  All  functions  of  X  and  M;l  are  tabulated  in  tables  I  and  II. 

Once  X  is  known  as  a  function  of  x,  the  following  functions  can 
be  found: 
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(22) 

(23) 


and 


-  ^  Q.  +  ~  Mi2  (-0.000112  -  0.0094k  +  0.4175)1  (33) 

e.^lSte 

tg  V  p 


5*  =  [a*  +  A  Mx2  (-0.0001X2  -  0.0094X  +  0.4175)J  (35) 


'Sr 


and,  for  assumed  values  of  0<  tj<-  1, 


u 

U1 


=  2n  - 


21)3  +  + 


=  1 


,  T-1  w  2 
+  ~o"  M1 


(36) 


(20) 


(30) 


The  velocity  and  temperature  profiles  can  then  he  found  in  terms  of  the 
physical  (x,y)  coordinate  system  hy  the  following  relation: 


* 
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where 


\2  o 

a2=k+lX+4 


X2 

a3  =  -  —  -  2X 


a  _  5Xf  .  4X 
a4  “  12  +  3  8 


5X2  5X 


a5  =  "~+"3+4 


a6  = 


5X' 

12 


3X  +  4 


a,  =  -  %  ♦  f  -  4 


_  X2  X  , 
a8  ~  36  "  3  +  1 


SPECIAL  CASES 


Solution  near  a  stagnation  point.  -  In  the  immediate  neighborhood 
of  a  stagnation  point,  uj_  is  equal  to  zero  and  the  incompressible 
solution  must  apply.  This  solution  is  (reference  3): 


XQ  =  7.052 


„  0.077 

0  duj_/dx 


(37) 


-0.0652 


d2u  jJ  dx^ 
dui/dx 


J 


Plat-plate  solution.  -  In  the  case  of  the  flat  plate,  X  =  0  a^d 
the  conditions  at  the  outer  edge  of  the  boundary  layer  are  constant  and 
equal  to  the  standard  condition(s) .  Equation  (27)  therefore  becomes 
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dZ 

dx 


=  ^  [fi°)  f(Ml)  +  f2(°0 


0.470 

u-, 


and 


Z  = 


0.470  x 

Ui 


From  equations  (25)  and  (34),  the  momentum  thickness  is 

Q  =  0.685 


1 


The  boundary- layer  thickness  is  found  from  equations  (22),  (34), 
and  (33)  to  be 

5  =  5.836 [l  +  0.417  Mx2  (^Jj 
The  displacement  thickness,  from  equations  (23)  and  (35),  is 
8*  .  1.75*/^  [l  +  1.392  (C±)  M^J 
and,  from  equation  (36),  the  vail  shearing  stress  is  found  to  be 


(38) 


(39) 


(40) 


(41) 


w  2.918  Vvix 


ul^l  x  /ulC 


(42) 


Velocity  and  temperature  profiles  can  be  found  as  functions  of  t]  by 
letting  \  vanish  in  equations  (20)  and  (30),  respectively.  The  rela¬ 
tion  of  rj  and  y  can  then  be  obtained  from  equation  (32).  The  error 
introduced  in  the  flat-plate  solution  by  assuming  that  the  velocity 
profile  is  a  fourth-degree  polynomial  can  be  deduced  by  comparing  these 
results  vith  the  results  of  reference  5: 
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Equations  (43)  and  (45)  apply  for  any  Prandtl  number,  whereas  equa¬ 
tion  (44)  applies  for  Pr  =  1  only. 

It  can  therefore  he  seen  that  the  results  of  the  present  analysis 
are  quite  accurate  for  the  case  of  zero  pressure  gradient.  No  definite 
conclusions  for  the  nonzero  pressure  gradient  case  can  he  made,  however. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  August  8,  1951 
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APPENDIX  A 


TRANSFORMATION  OF  MANGLER 


The  laminar-houndary- layer  equations  for  three-dimensional  flows 
with  axial  symmetry  are 


c)  /_  _ 
&  (pr°" 


pc£  pSy^Sy/ 

(Al) 

Jj  (pro7>  -  0 

(A2) 

<§>  +  5  /t:  StV  -  /  SuT 

S  3?  V  IwJ 

(A3) 

where  the  bar  is  used  to  differentiate  the  three-dimensional  from  the 
two-dimensional  quantities,  and  rQ(x)  defines  the  radial  coordinate  of 

the  body  in  a  meridional  plane.  The  continuity  equation  is  satisfied 
by  the  following  stream  function: 


r0 


(A4) 


r0  P  v  = 


The  following  transformation  variables  were  introduced  by  Mangier 
(unavailable  report): 


x  «  K£ 


0 


>2(x)  dx 


y  =  K  r0(x)  y 


(A5) 


where 


K  is  an  arbitrary  constant.  Thus 
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3  K2r  2  s  yV  5 

0  ^x  Tq 


where  the  prime  indicates  differentiation  with  respect  to  x. 

It  is  further  assumed  that 

p(x,y)  =  p(x,y) 
t(x,y)  =  t(x,y) 

P(x,y)  =  n(x,y) 
p(x,y)  -  p(x,y) 

?(x,y)  =  |  t(x,y) 

Equations  (Al) ,  (A2) ,  and  (A3)  can  then  he  transformed  to  the  following 
form: 


du 
U  chc 

,  du  1  dp 

+  =  "  p  ta 

+  - 
P 

a  /  au\ 

(A6) 

4  M  +  ^  (pv)  = 

0 

(A7) 

(  St 

pCP  (U  35  + 

at\  aP 
v^j  =  uS  + 

If 

1  dt\  /du\^ 

(A8) 

!? 

ll 

X 

2 

^y) 

r- 

f  "1 

(A9) 

v(x,y) 

-  55  [*<*,7)  - 

yro 

‘  Kr0 

2  u(x,y)| 

Equations  (A6),  (A7),  and  (A8)  are  now  identical  to  equations  (3),  (4), 
and  (5).  The  solution  of  these  equations,  as  presented  in  the  text  of 
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this  report,  therefore  applies  to  the  three-dimensional  problem,  hut 
the  final  results  must  he  retransformed  to  the  three-dimensional  quan¬ 
tities.  The  boundary- layer  thickness  8  is  not  defined  in  a  solution 
of  the  differential  houndary-layer  equations,  hut  its  definition  is 
required  for  the  solution  of  the  integral  equation.  The  transformation 
of  6  can  best  he  accomplished  by  examining  the  velocity  profiles: 

iL  =  JL  (l\ 

ui  ui  w 


j*.  =  3-(z\ 

ui 

From  equation  (A9),  however, 


u  u 


It  therefore  follows  that 


I  -  I 
S  B 

whence,  from  equation  (A5), 


8  =  Kr0(x)  8  (AlO) 


The  momentum  thickness  is  transformed  as  follows: 


(All) 
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and  similarly. 


&*(*)  =  5*(x) 

Tw(x)  =  Kr0  Tv(x) 


(A12) 


where  the  quantities  without  the  bar  (representing  two-dimensional 
quantities)  are  given  by  equations  (33)  to  (36). 


Therefore,  if  it  is  desired  to  determine  any  of  the  boundary- layer 
characteristics  at  a  distance  x  from  the  nose  of  an  axially  symmetric 
body,  the  corresponding  two-dimensional  distance  x  is  first  found  from 
equation  (A5).  The  two-dimensional  quantities  (&,  &*,  6,  «r  ,  and  so 

forth)  are  then  found  using  equations  (33)  to  (36),  and  are  transformed 
to  the  three-dimensional  quantities  at  x  through  the  use  of  equa¬ 
tions  (A10)  to  (A12).  The  transformation  applies  for  flows  with  or 
without  streamwise  pressure  gradients. 
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APPENDIX  B 


SOLUTION  OF  DIFFERENTIAL  EQUATION  BY  RUNGE-KUTTA  METHOD 

The  Runge-Kutta  method,  of  finite  differences  (reference  6)  has 
been  found  to  he  very  satisfactory  for  solving  equation  (27).  It  is 
believed  that  a  detailed  discussion  of  the  method,  as  applied  to  the 
present  problem,  will  be  of  help  to  the  reader. 


It  is  suggested  that  a  table  of  the  following  form  be  used: 


The  distance  along  the  body  (x)  is  subdivided  into  a  number  of  incre¬ 
ments  h,  and  the  free-stream  variables  u^,  du^/dx,  and  M^  at  each 

point  are  tabulated  in  columns  (2),  (3),  and  (4).  Column  (5)  follows 
from  equation  (28),  and  column  (6)  is  easily  tabulated.  Columns  (l) 
to  (6)  can  be  tabulated  for  all  x,  whereas  for  columns  (7)  to  (ll), 
each  row  must  be  calculated  separately. 

The  value  of  ZQ  =  ®o2/vs  (column  (7),  row  l)  is  given  by  equa¬ 
tion  (37)  for  stagnation-point  flow,  or  is  equal  to  zero  if  ©q  =  0. 
The  value  of  f^(X)  (column  8)  is  obtained  from  equation  (29),  or  in 
column  notation. 


fx(x)  =  (7) (3) (6) 
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Columns  (9)  and  (10)  are  obtained  from  figure  (l).  Column  (ll) 
lists  the  Runge-Kutta  parameter  k,  which  in  column  notation  is 


Rows  1  to  4  can  now  he  completed  in  this  manner.  At  the  completion  of 
row  4  the  total  increment  in  Z  can  he  computed  using  the  following 
equation: 


AZ  -  g  (kj_  +  2kg  +  2kj  +  k^) 

The  entire  process  is  then  repeated  for  rows  5  to  8.  The  value  of  Z ^ 

(column  (7),  row  5)  is  equal  to  Zq  +  AZ.  The  procedure  is  carried  out 

for  all  increments  of  x,  and  as  a  result  X  is  obtained  as  a  function 
of  x. 
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TABLE  II  -  FUNCTIONS  OF  M 


Mi 

ui 

1  +  Zjr  %2 

f(Mx) 

M1 

III  M  2 
2  M1 

1  +  ^  Mj2 

f(Mi) 

0 

0 

1.000 

-2.000 

2.2 

0.968 

1.968 

0.213 

.2 

.008 

1.008 

-1.974 

2.4 

1.152 

2.152 

.409 

.4 

.032 

1.032 

-1.861 

2.6 

1.352 

2.352 

.589 

.6 

.072 

1.072 

-1.698 

2.8 

1.568 

2.568 

.748 

.8 

.128 

1.128 

-1.489 

3.0 

1.800 

2.800 

.893 

1.0 

.200 

1.200 

-1.250 

3.2 

2.048 

3.048 

1.024 

1.2 

.288 

1.288 

-.994 

3.4 

2.312 

3.312 

1.141 

1.4 

.392 

1.392 

-.733 

3.6 

2.592 

3.592 

1.247 

1.6 

.512 

1.512 

-.476 

3.8 

2.880 

3.880 

1.343 

1.8 

2.0 

.648 

.800 

1.648 

1.800 

-.231 
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4.0 

3.200 

4.200 

1.500 
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Figure  2.  -  Concluded.  Auxiliary  functions  used  in  equations  (31)  and  (32). 
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bodies  with  axial  symmetry  through  the  use  of  bodies  with  axial  symmetry  through  the  use  of 

Mangier’ s  transformation.  Mangler's  transformation. 


